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et article vient en annexe des articles [B 9 734] et [BE 9 735] sur la theorie 
des systemes thermofrigorifiques. 

Le systeme thermofrigorifique etudie ici est un systeme a absorption continue 
utilisant le couple ammoniac-eau. 

II est destine a refroidir un flux de liquide frigoporteur pour le refroidissement 
indirect d'enceintes dont la temperature interieure est voisine de 0 °C. 

Nous allons d'abord presenter I'etude theorique d'un tel systeme puis nous en 
donnerons une application numerique. 

Nous utilisons comme unite de pression le bar absolu. Rappelons que : 

7 bar = 10 5 Pa 
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MACHINES THERMOFRIGORIFIQUES 


1.1 Calculs thermodynamiques 


Pour effectuer les calculs necessaires, il faut disposer de donnees. 

1 - Donnees imposees par la nature du probleme frigorifique a 
resoudre : raison d'etre, type de machine, implantation, etc. ; elles 
sont incontournables et I'ingenieur n'a aucun choix possible. Ces 
donnees sont : 

— la puissance frigorifique a produire qui peut etre fournie direc- 
tement, <1> F , ou etre deduite d'un bilan frigorifique effectue comme 
il a deja ete explique dans Particle de cette collection concernant la 
theorie des machines frigorifiques a compression [B 9 731, 
§2.4.1.1.41; 

— la temperature du milieu a refroidir : 

• temperature de I'enceinte froide (chambre, compartiment d'un 

meuble frigorifique, etc.), 

• temperature de sortie d'un fluide frigoporteur refroidi par le 

systeme frigorifique ; 

— la temperature d'entree du fluide de refroidissement dans les 
composants du systeme auxquels il faut extraire de la chaleur : 

• temperature du fluide de refroidissement du condenseur, 

• temperature du fluide de refroidissement de I'absorbeur. 

Bien que cela ne soit pas obligatoire, ce sont souvent les memes 
temperatures, le fluide de refroidissement etant I'eau ou plus rare- 
ment I'air (ce choix dependant, on s'en doute, de la puissance ther- 
mique a extraire). Ces temperatures dependent du lieu 
d'implantation de la machine. 

Notons que la temperature de la source de chaleur necessaire au 
chauffage du bouilleur peut etre, selon les cas, une donnee : 

— imposee, si la chaleur est recuperee d'une source ou la tempe- 
rature est connue (source de recuperation de chaleur issue d'un pro- 
cessus industriel donne...) ou limitee (energie thermique issue d'un 
capteur solaire, la temperature disponible etant tres variable selon 
le type de capteur...) ; 

— a se fixer, si I'energie thermique est non recuperee mais pro- 
duce pour les besoins du systeme. Une optimalisation peut alors 
etre necessaire. 

2 - Donnees que I'on doit se fixer ou estimer, le choix restant pos- 
sible puisque ces parametres ne sont pas fondamentalement impo- 
ses. Ces donnees sont : 

— la nature du couple frigorigene-absorbant ; en fait ce choix est, 
dans la pratique, extremement limite, comme nous I'avons vu en 
[BE 9 735, § 2.2] et, par ailleurs, selon la nature du probleme il se 
trouve, en fait, quasiment impose : pour le froid climatique, c'est 
pratiquement le systeme eau-bromure de lithium qui est choisi tan- 
dis que la production de froid pour atteindre des temperatures infe- 
rieures a 0 °C requiert I'emploi du couple ammoniac-eau ; 

— le titre en frigorigene de la vapeur issue du generateur lorsque 
la rectification s'impose ; 

— les « pincements » ou ecarts minimaux de temperature entre 
les fluides qui echangent de la chaleur dans les divers echangeurs 
thermiques de la machine. Le choix de ces « pincements » peut faire 
I'objet d'etudes d'optimalisation ; 

— les temperatures des points caracteristiques : 

• du circuit de frigorigene, 

• du circuit des solutions. 


Exemple : pour fixer les idees, considerons un systeme frigorifique 
simple a absorption continue du type de celui de [BE 9 735, figure 2] 
represente schematiquement par la figure 1 et destine au refroidisse- 
ment d'un debit de liquide frigoporteur. Celui-ci circule dans I'evapora- 
teur EV ou il se refroidit de d e F p a d sF p. Le condenseur CD du 
systeme est refroidi par une circulation d'eau qui s'echauffe de 0 e CD a 
6 S cd ■ De meme, I'absorbeur AB est egalement refroidi par de I'eau 
dont la temperature varie de 6 e AB a 6 sAB . Pour §tre complet, on sup- 
pose aussi qu'il s'agit d'un systeme ammoniac-eau, lequel impose une 
rectification. 

Le generateur de vapeur d'ammoniac, dont le titre en ammoniac est 
xpp~ 1 (a quelques milliemes pres), comprend done le bouilleur BO, la 
colonne de rectification CR et le dephlegmateur DE. On peut ainsi 
appliquer au frigorigene qui circule dans le circuit frigorifique CF les 
grandeurs thermodynamiques de I'ammoniac pur. 



Grandeurs 

p (bar) pression 
0(°C) temperature 
cji(K) temperature absolue 
x titre 

(kW) puissance thermique 
P SR (kW) puissance fournie 
a la pompe de solution riche 
m debit massique 

Indices 

C condensation ou chaud 
FR frigorigene 
e entree 
s sortie 
F froid 
r riche 
p pauvre 
FP frigoporteur 


Circuits 

CD condenseur 
BL bouteille de liquide 
DT detendeur 
EV evaporateur 
AB absorbeur 
PS pompe 
RS regleur 

ET echangeur thermique 
SP solution pauvre 
GE generateur de vapeur 
d'ammoniac 
CF circuit frigorifique 
CS circuit des solutions 
BO bouilleur 
SR solution riche 
CR colonne de rectification 
DE dephlegmateur 


Figure 1 - Schema d'un systeme frigorifique simple a absorption 
continue 
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MACHINES THERMOFRIGORIFIQUES 


1.1.1 Donnees necessaires 

Donnees imposees, sans choix possible 

Puissance frigorifique & F \ en donnant le debit massique m B p du 
liquide frigoporteur, de capacite thermique massique Cpp, a refroidir, 
on a alors la puissance frigorifique necessaire : 

0 F = mppc FP (0 e , FP - 9 S FP ) (W ou plus souvent kW) (1) 

Cette valeur, nette, doit etre, en general, legerement augmentee 
pour tenir compte des pertes thermiques du circuit frigorifique. 

Temperature de sortie 9 S FP du frigoporteur : c'est en effet une 
donnee essentielle du probleme frigorifique pose car c'est el le qui 
est demandee par I'utilisateur de la machine. 

Conditions de refroidissement du condenseur CD : il faut en effet 
connaTtre la temperature 9 e CD de I'eau entrant dans cet organe. 

Conditions de refroidissement de I'absorbeur AB : le probleme est 
identique a celui du condenseur ; il faut imperativement connaTtre la 
temperature d'entree d'eau 9 e AB dans cet echangeur ther-mique. 

Conditions de chauffage du bouilleur BO (qui peut etre impose ou 
non seion les cas) : si celui-ci est chauffe par une vapeur qui se 
condense (vapeur d'eau), on doit connaTtre la temperature de 
condensation 9Q V de la vapeur ou la pression de saturation p Bv cor- 
respondante. Si le bouilleur est chauffe par un fluide caloporteur 
(eau « surpressee », par exemple), on doit connaTtre la temperature 
de sortie 9 S b q de ce fluide. 

Donnees a fixer, done avec quelques possibilites de choix 

Titre x FR de la vapeur d'ammoniac produite : voisin de 1, il est 

necessaire a I'etude de la colonne de rectification. 

Pincement de I'evaporateur : c'est alors A0 p)n (9 S F p- 9 F ) ou 9 F 
est la temperature de vaporisation du frigorigene. Ce pincement 
etant fixe, la temperature 9 F se trouve determinee. 

Pincement du condenseur : c'est alors A 0 p ; a co = $C~ ®s,CDi ou 
9 s ,cd est l a temperature de sortie d'eau du condenseur et 9 C la tem- 
perature de condensation du frigorigene. Mais si I'on connaTt 9 e qd 
on ne connaTt pas, « a priori » coQtii depend de I'echauffement de 
I'eau dans le condenseur A9^ cR qd = (0 s ,CD _ ®e,CDi- On doit done 
aussi connaTtre ou se fixer cet echauffement. Lorsque A 9 p m qd et 
A Sech.CD sont fixes, alors la temperature de condensation 9q est 
connue. 

Pincement de I'absorbeur : c'est alors A 9 pmAB = ( 9 a - 9 sAB ) ou 9 a 
est la temperature de sortie de la solution riche de I'absorbeur et 
9 s ,ab temperature de sortie de I'eau de cet organe. Mais pour 
avoir acces a 9 S AB , on doit aussi, comme pour le condenseur, se 
fixer I'echauffement de I'eau de refroidissement 
Ae ech,AB = ( Qs.ab - e e,AB) ■ Lorsque A 9pin :AB et A 9ech,AB sont fixes, 
alors la temperature 9 a de sortie de solution riche est connue. 

Pincement du bouilleur : si celui-ci est chauffe par un fluide calo- 
porteur, comme le suggere la figure 1 , le pincement est alors 
A ®pin,BO= $s,BO-6d> ou 9 sB q est la temperature de sortie du fluide 
chauffant et 9 B la temperature de sortie de la solution pauvre. Pour 
connaTtre cette derniere, il est necessaire de se fixer aussi le refroi- 
dissement du fluide chauffant dans le bouilleur 

A9 ref bo = ( BO ~ bo ) ■ Si le bouilleur est chauffe par la conden- 
sation de vapeur, il sutfit, pour connaTtre 9 d , de se fixer le seul pince- 
ment qui est alors A9p in b q = ( 9cv~9d )■ 

Enfin, on doit aussi connaTtre, ou du moins estimer : 

— sur le circuit de fluide frigorigene, CF\ 

• la temperature 9 3 de sortie du liquide frigorigene de I'ensem- 

ble condenseur CD-bouteille de liquide BL. Cela revient a se fixer 

le sous-refroidissement sr= ( 9 C - 9 3 ) du liquide frigorigene, 

• la temperature 64 d'entree au detendeur DT si on pense que 

celle-ci est differente de la precedente, 



• la temperature 9 6 de la vapeur quittant I'evaporateur. Cela 
revient a se fixer la surchauffe s= (9g- 9 F ) de la vapeur en ce 
point, 

• la temperature de la vapeur de frigorigene entrant dans 

I'absorbeur. Cela revient a estimer I'echauffement 
A9 6c h v = cette vapeur dans la conduite d'alimenta- 

tion de I'absorbeur ; 

— sur le circuit des solutions, CS : 

• la temperature 9 a etant determinee comme on I'a indique plus 
avant, il nous faut estimer les imperfections de I'absorbeur pour 
connaTtre le point a' du cycle « reel » des solutions conforme- 
ment a la figure 2 , done le titre reel en frigorigene x r de la solution 
riche (ici x' r ), 

• la temperature 9g de sortie de la pompe PS est evidemment un 
peu superieure a la temperature 9 a en raison de la puissance p PSr 
fournie a la pompe et des imperfections de celle-ci. Comme on ne 
connaTt pas Pp Bn on ne peut connaTtre immediatement 9 B . 
Neanmoins, compte tenu du debit de solution riche, la difference 
de temperature (9 B - 9 a ) est toujours petite et peut etre, de prime 
abord, ignoree, 

• la temperature 9 C de la solution riche entrant dans la colonne 
de rectification (designee 9 ( SR dans [BE 9 735, figure 18]) est 
prise, des I'abord, comme etant celle du point de bulle de cette 
solution sous la pression de la colonne p B (ou pq + 8pq b , confor- 
mement a la figure 2 , s'il doit exister une chute de pression sensi- 
ble entre le generateur GE et le condenseur CD - point b' de cette 
figure), 

• les temperatures 9 e et 9f ne peuvent etre « a priori connues ». 
Elies devront etre determinees ulterieurement par les calculs ther- 
modynamiques afin notamment de verifier que 9 e > et que 
I'ecart entre ces temperatures est suffisant, 

• le cas echeant, on estimera I'imperfection du generateur 
(bouilleur) (figure 2 ) pour determiner le titre reel en frigorigene x p 
de la solution pauvre (ici x' p ). 

1.1.2 Marche a suivre 

On effectue successivement les operations suivantes : 

— determination du titre x r de la solution riche, apres evaluation 
de la temperature 9 a et estimation de I'imperfection de I'absorbeur ; 
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x 


Figure 3 - Cycle des solutions dans le diagramme 
enthalpie-concentration 


— determination du titre x p de la solution pauvre, apres evalua- 
tion de la temperature 8 d et estimation de I'imperfection du gene- 
rateur (bouilleur) ; 

— etablissement du tableau des grandeurs thermodynamiques 
du frigorigene (ammoniac) pour les differents points de son cycle : 
pression, temperature, enthalpie massique, volume massique. Les 
pressions et temperatures ayant ete estimees comme indique ci- 
avant, ce tableau s'etablit sans difficulty ; 

— debut du trace du cycle des solutions dans le diagramme 
enthalpie-concentration h-x : 

• on trace I'isotitrique x r en portant sur celle-ci (figure 3 ) le point 
a, defini par sa pression et sa temperature, et le point c, defini de 
la meme maniere, 

• il n'est pas encore possible de placer le point b, I'energie 
requise par la pompe n'etant pas encore determinee, ce qui ne 
permet pas de connaitre I'echauffement (faible) de la solution, 

• on trace I'isotitrique x p en portant sur celle-ci le point d, 

• il n'est pas encore possible de placer les points e et f. II faut, 
pour cela, attendre d'avoir fait le bilan energetique de I'echangeur 
des solutions El ; 

— etablissement du tableau des grandeurs thermodynamiques 
des solutions pour les divers points de leur cycle : titre, temperature, 
enthalpie massique. Ce tableau ne peut etre encore qu'imparfaite- 
ment rempli puisque des donnees manquent encore ; 

— execution des calculs thermodynamiques concernant les 
divers composants. On tire de ces calculs les grandeurs thermody- 
namiques encore inconnues, ce qui permet de completer les 
tableaux et le cycle des solutions ; 

— execution de I'epure de la rectification et determination finale 
des puissances thermiques <P M et <P DE faites apres I'etude de la 
colonne de rectification (§ 1.2) ; 

— determination du coefficient de performance du systeme et 
comparaison avec celui du systeme tritherme ideal ; 

— determination des diverses conduites du systeme ; 

— determination des debits des fluides de refroidissement ; 

— determination du debit de fluide de chauffage du bouilleur. 


1.1.3 Relations thermodynamiques 

L'application du premier principe de la thermodynamique a un 
« systeme ouvert » (or tous les composants du circuit frigorifique : 
evaporateur, condenseur, absorbeur, pompe etc., pris individuelle- 


ment, sont des « systemes ouverts » avec entree et sortie de flui- 
des), donne, par unite de masse de fluide en circulation : 


w+ q = (h s - h e ) + \AE p \ s e +\A E c \ s g + ... (kJ/kg) (2) 


avec w travail mecanique echange entre le systeme et 

I'exterieur, 

q energie thermique echangee entre le systeme et 

I'exterieur, 


h enthalpie massique (e, pour I'entree et s, pour la 

sortie du systeme), 

AE p |® variation d'energie potentielle gravifique du 
fluide entre e et s (energie exterieure), 

AE C |^ variation d'energie cinetique du fluide entre e et 
s (energie exterieure). 

Nota : on se limitera a ces deux formes d'energie exterieure ; mais il pourrait y en avoir 
d'autres. 


Comme c'est habituel dans nos types de systemes frigorifiques, 
on peut, tout a fait, negliger les energies potentielle gravifique et 
cinetique. La relation (2) se simplifiant alors : 


w+ q = (h s - h e ) (kJ/kg) 


( 3 ) 


On applique cette relation aux divers composants du systeme. 

Debit massique de frigorigene m F en circulation 

Dans I'evaporateur EV, comme w= 0 (echangeur statique), on a, la 
detente du frigorigene dans DT etant isenthalpique : 



& F = m F (h 6 -h s ) = m F (h 6 -h 4 ) 

(kW) 

( 4 ) 


&F 



d'ou : 

m F = — 

(kg/s) 

( 5 ) 


(h 6 -h 4 ) 




Puissance thermique & CD echangee au condenseur 

Le meme debit m F traverse, en regime permanent, le condenseur 
CD. D'apres (3), la puissance thermique echangee dans cet organe 
est alors : 


0> c = m F (h 2 - h 3 ) (kW) (6) 

Debits volumiques de frigorigene 

Debit volumique de frigorigene vapeur V F v, 7 entrant, en 1, dans 
I'absorbeur : 


Vfv. i = Fb F v FV , (m 3 /s) ( 7 ) 

avec v F v i volume massique de la vapeur de frigorigene en 1. 

Debit volumique de frigorigene vapeur Vfv, 2 sortant, en 2, du 
generateur : 

Vfv, 2 = mFV FV ,2 (m 3 /s) (8) 

avec v F v ,2 volume massique de la vapeur de frigorigene en 2. 

Debit volumique de liquide frigorigene Vfl ,4 entrant, en 4, dans 
le detendeur DT : 


Vfl .4 = m F v FL4 (m 3 /s) ( 9 ) 

avec Vp V ,4 volume massique du liquide frigorigene en 4. 

Taux de circulation et debits des solutions 

— taux de circulation de la solution riche : en reprenant les bilans 
exposes en [BE 9 735, § 2.5], mais en tenant compte du fait que la 
vapeur de frigorigene qui circule dans CE n'est pas pure mais a un 
titre xppj < 1 (tout en etant cependant tres voisin de 1), on a : 

x fr ~ x p 

T c r = kg/s de solution riche par kg/s de frigorigene ( 10 ) 

X r~ x p 
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— debit massique msR de solution riche : 

msR = f c r mp (kg/s de solution riche) (11) 

— debit volumique Vsr jS de solution riche sortant de I'absor- 
beur en a : 

VsR,a= m S RV S R a (m 3 /s de solution riche) (12) 

II faut done disposer de donnees concernant les volumes massi- 
ques (ou les masses volumiques) des solutions frigorigene-absor- 
bant, pour diverses concentrations en frigorigene et dans le 
domaine des temperatures rencontrees. 

— taux de circulation de la solution pauvre : d'apres [BE 9 735 , 
relation (8)], on a : 

X FR ~ x r 

T = (kg/s de solution pauvre par kg/s de frigorigene) (13) 

• R x r -x p 

— debit massique msp de solution pauvre : 

msp = i c p mp (kg/s de solution pauvre) (14) 

— debit volumique Vsp, d de solution pauvre sortant du bouilleur 
en d : 

Vsp ,d = mspVsp d (m 3 /s de solution pauvre) (15) 

Puissance de la pompe PS de solution riche 

Cette pompe deplace le debit volumique de solution riche. I/s/?, a , 
de I'absorbeur, ou regne la pression p F , vers le generateur oil la 
pression est p E . Si son fonctionnement etait ideal, la puissance 
mecanique qu'elle consommerait serait : 

P ps, i = VsR r a{Pc~ P f) (W) (16) 

Mais on doit tenir compte de I'imperfection du groupe moto- 
pompe en admettant, pour celui-ci, un rendement energetique rai- 
sonnable r] MP . On a alors la puissance reelle du moteur de cette 
pompe : 

Vsr, a(Pc~ Pf) 

Pp S ,r= Jf (W) (17) 


ou Cg R est la capacite calorifique massique de la solution riche a la 
sortie de I’absorbeur. Comme Ah est faible, la variation de tempera- 
ture AS est des plus reduites. 

Puissance thermique echangee dans I'echangeur thermique El 
intersolutions 

Toujours d'apres la relation (3), en admettant que I'echangeur El 
soit parfaitement isole thermiquement de I'exterieur : 

&EI = m SR<h c - h b> (kW) (20) 

Refroidissement de la solution pauvre dans I'echangeur El 
intersolutions 

Si I'on considere comme nuls les echanges thermiques de cet 
echangeur El avec I'exterieur (en realite ces echanges sont faibles 
car El est isole thermiquement) on a, en appliquant la relation (3) : 

‘Pei = rnsp(h d - h e ) (21) 

®EI 

d'ou : h g = h d ~— — (kJ/ kg de solution pauvre) (22) 

msp 

On en deduit la temperature 8 e du point e et la faisabilite de 
I'echangeur El, faisabilite qui est assuree si 8 e est suffisamment dif- 
ferente de 9 b . On peut maintenant poursuivre le trace du cycle des 
solutions (figure 3) jusqu'en ce point, evidemment situe sur I'isoti- 
trique x p . 

Determination de I'enthalpie h f a la sortie du regleur RS 

Appliquons la relation (3) au regleur de solution pauvre RS. Pour 
ce composant w= 0 (e'est un robinet, parfaitement statique) et q = 0 
(il n'y a pas d'echange thermique, la surface d'echange thermique 
offerte etant negligeable). De ce fait, (hf - h e ) = 0 , la detente de la 
solution est isenthalpique et hf= h e . Sur le cycle de la figure 3, les 
points e et f sont confondus. 

Bilan energetique de I'absorbeur. Calcul de ® AB 

+ m F h 1 -® AB - m SR h a + m S ph f = 0 (23) 


Determination de I'enthalpie h b et de la temperature de la 
solution riche a la sortie de la pompe 

En admettant que les echanges thermiques entre la pompe et 
I'environnement soient negligeables, alors la variation d'enthalpie 
est, d'apres (3) : 


On peut done deduire, de ce bilan, la puissance thermique a 
extraire de I'absorbeur : 

<P A B = m F h 1 + msph f -m S Rh a (kW) (24) 


( h b~ h a) = w = 


PS. r 


m sr 


(kJ/kg de solution riche) 


variation d'enthalpie massique generalement tres faible. 


(18) 


Dans la relation (18), on a utilise la puissance reelle, ce qui 
sous-entend que les pertes de la motopompe sont integrale- 
ment absorbees par le fluide qui circule. Ce n'est pas absolu- 
ment exact. En effet, les pertes du moteur electrique se 
retrouvent dans I'air et non dans la solution en circulation. La 
valeur trouvee en (18) est done une valeur par exces. L'erreur 
que I'on commet ainsi est minime. 

Sur le trace du cycle des solutions, les points a et b sont pratique- 
ment confondus. La temperature 6 b est tres peu differente de S a car : 

h b~ h a 

(0 b -S a ) = — (K) (19) 

C SR 


Bilan energetique du generateur. Calcul de I'ensemble (O m - O de ) 

-m F h 2 +® M - <P DE - msph^F m SR h c = 0 (25) 

On peut done deduire, de ce bilan, la puissance thermique globale 
a fournir au generateur. Celle-ci est la somme algebrique de la puis- 
sance thermique a apporter au bouilleur, ® M , et de la puissance 
thermique a extraire (done negative) au dephlegmateur, ® DE : 

( 0 m~ 0 de) = m Fh 2 + m S ph cl -msRh c (kW) (26) 

Verification du premier principe de la thermodynamique 
applique a I'ensemble du systeme 

L'ensemble des calculs des diverses puissances mises en jeu dans 
les composants du systeme a absorption considere est desormais 
termine. On doit maintenant s'assurer que le premier principe de la 
thermodynamique applique a cet ensemble, pris comme etant un 
systeme ferme, est bien verifie. II faut, pour cela, que les puissances 
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regues de I'exterieur par I'ensemble du systeme et les puissances 
qu'il cede a I'exterieur soient egales : 

®F+ (®M“ ®Df) + Pps,r= ( t‘CD + ®AB (27) 

Nous ne pouvons pas aller plus loin, maintenant dans nos 
calculs. En effet, ayant affaire a un systeme ammoniac-eau, 
I’emploi d'une colonne de rectification s'impose afin d'epurer la 
vapeur ammoniacale de I'eau qu'elle renferme. C'est de I'etude 
de cette colonne que I'on tirera les valeurs thermodynamiques 
qui nous manquent encore : 

— la puissance thermique motrice & M a fournir au bouilleur 
BO; 

— la puissance thermique 0q E a extraire au dephlegmateur 
DE ; 

— le coefficient de performance frigorifique COP F e t le rende- 
ment de cycle T ] c ; 

— les debits de fluides chauffant et refroidissant ainsi que les 
conduites necessaires. 

Nous appuyant sur ce qui a ete expose dans [B 9 735, § 2.3 et 
2.7] nous appliquerons la methode graphique (mais rigoureuse) 
de Ponchon et Savarit a I’etude de la colonne de rectification. 


1.2 Etude d'une colonne de rectification. 
Application de la methode 
de Ponchon et Savarit 

Dans ce qui suit, nous supposerons qu'au niveau de chaque pla- 
teau de la colonne la vapeur qui le quitte est en equilibre thermody- 
namique (temperature et concentration en frigorigene) avec le 
liquide qui en sort (rendement d'echange de chaleur et de matiere 
egaux a I'unite, ce qui n'est qu'une approche, bien sur). 

Donnees necessaires 

On se fixe, ou la nature du probleme a resoudre fixe : 

— la pression p qui regne dans la colonne. C'est la pression de 
condensation ou du moins une pression un peu superieure du fait 
des pertes de charge. El le permet de selectionner, dans le dia- 
gramme enthalpie-concentration h-x, la courbe « de bulle », B, la 
courbe « de rosee », R et la courbe auxiliaire CA correspondante 
(figure 4 ) ; 

— le titre x FR en frigorigene de la vapeur V FR sortant de la 
colonne. II fait I'objet d'un choix arbitraire ou d'une optimalisation ; 

— le titre x r en frigorigene de la solution riche liquide qui entre 
dans la colonne. Ce titre depend, comme nous I'avons vu plus haut : 

• de la pression d'evaporation p F , 

• de la temperature de I'absorbeur, 

• des imperfections de I'absorption ; 

— le titre x p en frigorigene de la solution pauvre liquide qui sort 
du bouilleur. Ce titre depend : 

• de la pression dans la colonne, 

• de la temperature du bouilleur, 

• des imperfections de celui-ci. 

Mode operatoire 

Dans le diagramme enthalpie-concentration (figure 4 ), on trace 
les isotitriques : 

— x FR , ce qui donne le point V FR sur la courbe « de rosee » R, 
I'etat thermodynamique du liquide de reflux L Re en equilibre avec 
cette vapeur souhaitee et la temperature 0 e Re de ce liquide (mais 
0 f Re peut etre une donnee liee au refroidissement du dephlegma- 
teur. x FR peut alors en dependre) ; 

— x r , ce qui donne le point L ER sur la courbe « de bulle » B si I'on 
admet que ce liquide est saturant. On determine la temperature 
d'equilibre 8 e SR de ce liquide sous la pression de la colonne ; 


— x p , ce qui donne le point L EP sur la courbe « de bulle ». Par la 
construction classique utilisant la courbe auxiliaire CA, on dete- 
rmine I'isotherme 0 t sp (temperature du bouilleur) et I'etat thermo- 
dynamique l/ 7 de la vapeur quittant le bouilleur pour entrer dans la 
colonne de rectification, vapeur supposee en equilibre avec le 
liquide sortant L$ P . 

En prolongeant I'isotherme 0 e SR au-dela de la courbe « de 
rosee », jusqu'a ce qu'elle rencontre I'isotitrique x FR (droite tracee 
en pointille), on obtient le point P 7(m , position limite inferieure du 
pole de purification. De meme, en prolongeant cette isotherme dans 
I'autre sens, au-dela de la courbe « de bulle », jusqu'a son intersec- 
tion avec I'isotitrique x p , on obtient le point Pj im , position limite 
superieure du pole d'epuisement. La droite principale de rectifica- 
tion PP' doit avoir une pente superieure a celle de la droite 
P ij m . P'nm. ■ Elle doit egalement passer par le point caracteristique L ER 
correspondant a la solution riche entrante. Au point P 7jm corres- 
pond, selon le diagramme de la figure 4 , I'extraction de chaleur 
minimale egale a QoE.lim. (kJ/kg) de frigorigene. Pour que la pente 
de la droite PP' des poles ait une pente superieure a celle de I'iso- 
therme Op SR il faut extraire du dephlegmateur une quantite de cha- 
leur par unite de masse de frigorigene O de > OoE.lim.- On pose done 
O de = kO DE n m avec k> 1. Le facteur /rest un parametre d'optimisa- 
tion. On obtient ainsi le point P et on peut tracer la droite des poles 
PP' qui passe par le point L ER correspondant a la solution riche 
entrant dans le generateur (c'est le point c de la figure 1 ). 

Le point L$ P , correspondant a la solution pauvre quittant le 
bouilleur (c'est le point d de la figure 1), etant connu, la quantite de 
chaleur Q M a fournir au bouilleur par unite de masse de solution 
pauvre sortante, est representee par le segment P' L SP . 

Conformement a la theorie exposee ci-avant [B 9 735], en com- 
mengant, par exemple, le trace a partir du plateau inferieur de la 
colonne, on joint par une droite P' et V F L'intersection de P' V 1 et 
de la « courbe de bulle » donne le point i 7 representatif de I'etat du 
liquide quittant le 1 er plateau. La radiale P'L 1 V 1 correspond au 
plateau © de la colonne d'epuisement. Si, comme on I'a suppose, la 
vapeur V 2 quitte ce 1 er plateau en equilibre thermodynamique avec 
le liquide L 7 qui en sort, il suffit, pour obtenir V 2 de tracer, a partir de 
L 7 , I'isotherme 0 7 par la construction que nous connaissons ; I'inter- 
section de cette isotherme et de la « courbe de rosee » donne le 
point V 2 . Si le point i/ 2 est a gauche de la droite PP' , comme c'est le 
cas sur la figure 4 on trace la droite P'V 2 , ce qui donne, sur la 
courbe B, le point L 2 correspondant au liquide en provenance du 
2 e plateau. La radiale P'L 2 V 2 correspond au plateau @. La cons- 
truction, a partir de L 2 de la nouvelle isotherme 0 2 nous conduit au 
point V 3 qui est maintenant a droite de la droite PP' des poles. II 
nous faut done utiliser a present, pour continuer le trace, le pole de 
purification P. A cela pres, le trace se deroule de la meme fagon : 

— on joint Pet V 3 , ce qui donne la radiale PV 3 L 3 e t le point L 3 ; 
cette radiale correspond au plateau (3). On construit alors I'iso- 
therme 0 3 ce qui donne le point V 4 ; 

— on joint Pet V 4 , ce qui donne la radiale PV 4 L 4 et le point L 4 ; 
cette radiale correspond au plateau @. On construit maintenant 
I'isotherme 0 4 , ce qui donne le point Vg ; 

— on joint Pet V s , ce qui donne la radiale PV s L s e t le point L s ; 
cette radiale correspond au plateau (5). On construit I'isotherme 0g, 
ce qui donne le point l/gqui correspond a la phase gazeuse entrant 
dans le dephlegmateur. 

L'ultime radiale PV 6 donne le point L' Re qui correspond au liquide 
de reflux sortant du dephlegmateur DE. La phase gazeuse V Fs quit- 
tant ce dephlegmateur, et supposee en equilibre avec L' Re , est un 
peu plus riche en frigorigene que la donnee de base x Fs > x FR . 

La colonne comprend done, entre le bouilleur BO et le deplegma- 
teur DE, 5 plateaux theoriques (parfaits), 2 dans la partie epuise- 
ment et 3 dans la partie purification. Si I'on tenait compte de leur 
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V FR en fait 

V Fs > V FR 



x Fs 

Figure 4 - Application de la methode de Ponchon-Savarit dans le diagramme enthalpie-concentration 


Toute reproduction sans autorisation du Centre frangais (Sexploitation du droit de copie est strictement interdite. ~i 

©Techniques de I'lngenieur, traite Genie energetique BE 9 73G - 7 


Ce document a ete delivre pour le compte de 7200092269 - cerist // 193.194.76.5 


MACHINES THERMOFRIGORIFIQUES 


imperfection par des rendements de plateau, le nombre reel serait 
un peu superieur. 

Nota : la solution riche arrive sur le plateau © avec une 
teneur en frigorigene un peu plus grande que le liquide L 2 . C'est 
une petite cause d'irreversibilite thermodynamique. 


Quelques remarques pratiques 

1 - Dans la colonne d'epuisement, on a une succession de temperatu- 
res toutes plus elevees que celle, 9 e SR , de la solution pauvre entrante. 
En descendant le long de la colonne le liquide s'echauffe et cede, par 
vaporisation du frigorigene : le liquide « s'epuise » en frigorigene. 

2 - Dans la colonne de purification, on a une succession de temperatu- 
res toutes plus basses que 9 ( SR . En remontant cette colonne on rencon- 
tre un liquide de plus en plus froid qui condense d'avantage la substance 
absorbante : la vapeur « se purifie » en frigorigene. 

3- Plus on accroit Qq E (figure 4), c'est-a-dire le reflux, plus le point P 
s'eloigne de P lim.. plus la droite des poles, PP ' , se redresse et s'eloigne 
de I'isotherme 9 ( SR , plus les ecarts thermodynamiques au niveau des 
plateaux augmentent, plus la colonne se simplifie car le nombre de pla- 
teaux necessaires diminue. 

4 - Par contre, plus on accroit Qq E , plus on augmente O ^ la puissance 
thermique necessaire au bouilleur, ce qui diminue le coefficient de perfor- 
mance de la machine. 

5 - Si le liquide riche entrant dans la colonne en E n'est pas dans les 
conditions de saturation correspondant a la pression p, done au point de 
bulle 9 € SR , (point L ER de la figure 5 ) mais a une temperature inferieure 
(point L e i de la figure 5 correspondant a un liquide sous-refroidi), rien 
n'est fondamentalement change dans ce qui precede. La construction de 
la droite P 7 P 7 ' , qui passe par L E i s'effectue de la meme maniere et est 
soumise aux memes regies. Cela etant, on constate, a pentes egales de 
P 7 P/ et de PP' : 

— qu'il faut accroitre I'apport calorifique au bouilleur Qm,i> Qm> ce 
qui est defavorable ; 

— que I'extraction de chaleur au dephlegmateur est plus reduite : 
@DE, i < Qde > 

— que I'importance relative de la colonne d'epuisement (dont les 
radiales des plateaux sont comprises entre P\ V 1 et la droite des poles 
P7P7 ) augmente par rapport a celle de la colonne de purification (dont 
les radiales s'echelonnent entre P 7 P 7 et P-iL Re ). 

6-11 est tout a fait important de chercher a reduire les irreversibilites 
thermodynamiques en s'attachant a minimiser, autant que faire se peut, 
les ecarts de temperature entre les milieux qui echangent de la chaleur. 
Ainsi la figure 6 donne le schema d'un generateur et de sa colonne de 
rectification elabores pour appliquer ce principe : 

— dans I'echangeur de solution pauvre ESP, lie aux plateaux les plus 
bas de la colonne d'epuisement CE, la solution pauvre SP commence a se 
refroidir de 1 a 2 en remontant dans la colonne avant d'entrer dans 
I'echangeur thermique intersolutions El 011 elle se refroidit de 2 a 3. El le 
est ensuite introduite dans I'absorbeur par le regleur de solution, RS. Le 
liquide qui descend dans cette partie de la colonne est done prechauffe 
par ESP avant d'entrer dans le bouilleur BO ; 

— a I'autre extremite, le dephlegmateur DE est precede de I'echan- 
geur de solution riche ESR lie aux plateaux les plus hauts de la colonne de 
purification CP. La solution riche et froide y entre en a et s'y echauffe 
jusqu'en b tandis qu'elle contribue a refroidir la vapeur avant que celle-ci 
entre dans le dephlegmateur DE. La solution riche SR, prechauffee, entre 
dans I'echangeur El pour s'echauffer de b a c avant son entree dans la 
colonne de rectification. 

II va de soit qu'on doit optimaliser I'ensemble des trois echangeurs 
ESP, ESR et El. 

7 - On peut ecjalement appliquer a I'absorbeur les memes principes 
que Ton vient d'enoncer. Ainsi dans I'absorbeur a ruissellement schema- 
tise dans la figure 7, la solution pauvre s'ecoule d'abord sur I'echangeur 
de solution riche de I'absorbeur ESRA parcouru interieurement par la 
solution riche et froide puisee au fond de la soute de I'absorbeur AB et qui 
remonte vers la pompe de solution PS. En collectant une partie de la cha- 
leur d'absorption du frigorigene, la solution riche s'echauffe de 1 a 2. La 
solution pauvre qui s'enrichit ruisselle ensuite sur la surface, plus froide, 
du refroidisseur R parcouru par le fluide de refroidissement, generale- 
ment de I'eau (on note le sens d'ecoulement du liquide de refroidisse- 
ment, a contre-courant par rapport a I'ecoulement de la solution 
ruisselante). Ce nouvel echangeur ESRA doit aussi etre pris en compte 
dans I'optimalisation globale du circuit des solutions entre I'absorbeur, 
I'echangeur intersolutions et le generateur dans son ensemble. 



1.3 Calculs thermodynamiques 
complementaires 

L'epure de la colonne de rectification etant terminee, il est possi- 
ble, maintenant de determiner les grandeurs thermodynamiques 
qui nous manquaient. 


1.3.1 Determination des puissances thermiques 
et ®DE 

Determination de 

L'epure de la figure 4 permet de mesurer le segment P'L sp qui 
represente, a I'echelle de cette epure, la quantite de chaleur Q M (kJ/ 
kg de solution pauvre en circulation) a apporter au bouilleur. Le 
debit massique de cette solution pauvre etant msp, d'apres la 
relation (14), on a : 

0 M = Q M m S p (kW) (28) 

Determination de & DE 

De la meme fagon, la quantite de chaleur, Q DE (kJ/kg de frigori- 
gene en circulation) a extraire au dephlegmateur, est determinee, a 
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Vers le condenseur 



De I'absorbeur Vers I'absorbeur 


Figure 6 - Schema d'un generateur et de sa colonne de rectification 


Vers I’echangeur 
thermique 

De I’evaporateur intersolutions El 



Figure 7 - Schema d'un absorbeur a ruissellement 


1.3.3 Coefficient de performance du systeme 
tritherme ideal - Rendement de cycle 


II est tres utile de comparer la valeur ainsi trouvee du coefficient 
de performance a la valeur calculee pour un systeme tritherme ideal 
fonctionnant entre les memes temperatures. Celui-ci est, d'apres 
[BE 9 734, equation (1)] (T etant les temperatures absolues) : 


copp 


(id) 


T f (T m - T c ) 
(Tc~ T f ) T m 


Le rendement de cycle I 7 c du systeme reel est done : 
COP 


copp iid) 

I caracterise le degre de qualite du systeme reel. 


(31) 


1.3.4 Autres calculs 


I'echelle du diagramme, par le segment PV FR . Pour le debit m F de 
frigorigene deduit de la relation (5), on a : 

&de = 0 DE m F (kW) (29) 

On s'assure ensuite que la valeur de (4>(vt _ ®D£) deduite de ces 
relations verifie bien le resultat obtenu par I’equation (26). 


1.3.2 Coefficient de performance du systeme 

Le coefficient de performance est : 

Effet utile 

COP = = (30) 

Ce que Ton paie &M + P PS,r 

Notons que la puissance mecanique P p $ r que I’on apporte au sys- 
teme pour le fonctionnement de la pompe est extremement faible, 
souvent pratiquement negligeable, devant la puissance thermique 
& M fourr\\e au bouilleur. 


Calcul des conduites necessaires 

Les divers debits volumiques : 

— du frigorigene dans les etats liquide et vapeur ; 

— des solutions liquides frigorigene + absorbant ; 

ayant ete determines, il suffit de choisir une vitesse convenable pour 
ces divers fluides afin d'en deduire les diametres interieurs des 
conduites necessaires. 

Exemple : on peut choisir comme vitesse de fluide entre 8 et 15 m/s 
pour les phases « vapeur » et entre 0,5 et 2 m/s pour les phases 
« liquides ». Pour le frigorigene liquide, entre les points 3 et 4, il 
convient d'adopter une vitesse de circulation particulierement reduite 
pour eviter la formation de la phase gazeuse, par chute de pression, au 
sein de ce liquide ordinairement tres peu sous-refroidi. 

Calcul des debits volumiques des fluides de refroidissement du 
condenseur, de I'absorbeur et du dephlegmateur 

On a la forme generale suivante pour le composant / : 

= V eau , j p ggu c eau (6 s j- 0 e j) (kW) (32) 
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et : V eau j= — (m 3 /s) (33) 

Peau^eau^s, i~ ®e, /) 

Calcul du debit de fluide chauffant du bouilleur 

S'il s'agit d'un liquide chauffant, son debit est deduit de relations 
analogues aux dernieres. Si le chauffage du bouilleur se fait par 
condensation de vapeur, le debit de celle-ci est calcule a partir de 
son enthalpie de condensation L VE : 

mvE = (kg/s de vapeur) (34) 

l ve 


On se reportera a la figure 1. 


2.1 Donnees 

Conditions de fonctionnement de I'evaporateur 

II doit refroidir de 8 e F p = - 5 °C a d s ,pp = - 10 °C un debit volumique 

de Vpp = 30 m 3 /h d'un liquide frigoporteur constitue par un 
melange eau + ethanol dont le point de congelation commengante 
est de - 20 °C. 

On admet un « pincement » de I'evaporateur de 5 K. La tempera- 
ture d'evaporation de I'ammoniac est done 8p = -15 "C. A cette tem- 
perature correspond la pression de saturation de I'ammoniac de 
Pl t= 2,36 bar [1]. Le detendeur d'alimentation de I'evaporateur en 
frigorigene DT permet de maintenir a la sortie de I'evaporateur, en 
6, une surchauffe de s = 8g- 9 f =4K. La temperature du point 6 est 
done 8g = - 11 °C. 

Puissance frigorifique a produire 

Le melange liquide frigoporteur renferme 29,5 % en masse 
d'ethanol. Sa masse volumique moyenne (a - 8 °C) est p F p = 969 kg/ 
m 3 et sa capacite thermique massique sous pression constante est 
c FP = 4 167 J/(kg ■ K) [2]. La puissance frigorifique nette a produire 
est done : 

< t>F = ^FPPppCpp{8 g pp- 9 s ,pp) (35) 

30 4 167 

& F = x 969 x x 5 = 1 68,2 kW 

3 600 1 000 

Majorons ce resultat de 5 % pour tenir compte des pertes thermi- 
ques du circuit frigorifique soit : 

0 F = 168,2 x 1,05= 176,7 

disons : 

0 F = 177 kW 

Titre massique en ammoniac du frigorigene 

II devra etre egal, au moins, a x F p = 99,5 % 

Conditions de fonctionnement du condenseur 

La temperature de I'eau disponible est 9 e CD = 20 ’C. Admettons 
un echauffement de I'eau dans le condenseur de A8 actl q d = 5 K. La 


temperature de sortie d'eau est 9 S q d = 25 °C. Fixons-nous un pince- 
ment du condenseur de A 9 pm cd= 6 K. La temperature de condensa- 
tion de I'ammoniac est alors 8q=30°C. A cette temperature 
correspond la pression de saturation de I'ammoniac de 
Pc= 11,67 bar [1], La temperature du liquide frigorigene sortant, en 
3, du condenseur est comprise entre la temperature de condensa- 
tion et la temperature d'entree d'eau au condenseur. Admettons que 
le sous-refroidissement du liquide forme soit sr= (9 q - 83) = 4 K. On 
a alors 83 = 26 °C. Cette temperature etant voisine de la temperature 
ambiante, nous admettrons que dans la conduite alimentant le 
detendeur DT la temperature du frigorigene ne varie pas 
84=83= 26 “C. 

Conditions de fonctionnement de I'absorbeur 

On utilise la meme eau et 8 e AB = 20 °C. Pour ameliorer I'absorp- 
tion, on s'attache a refroidir le plus possible I'absorbeur et on admet 
done un echauffement d'eau plus reduit Ad^ A b = 4 K ainsi qu'un 
pincement egalement plus faible A 8 p j n A b = 4 K. La temperature de 
sortie de solution riche de I'absorbeur est alors 8 a = 28 °C. Par 
ailleurs, compte tenu de la longueur et de I'isolation thermique de la 
conduite reliant I'evaporateur a I'absorbeur, on estime que la vapeur 
s'y rechauffe de A f?ec/i,v= (#j - 8 g) = 6 K. La temperature de la vapeur 
de frigorigene entrant dans I'absorbeur est done = - 5 °C. Pour 
simplifier, on supposera que la chute de pression entre evaporateur 
et absorbeur est negligeable. Ce dernier est done sous la pression 
p F = 2,36 bar. 

Conditions de fonctionnement du bouilleur (generateur) 

II est alimente en puissance thermique par de la vapeur d'eau 
sous la pression absolue de 3 bar (temperature de saturation 
=133 °C) qui provient de la chaufferie de I'usine oil les enceintes froi- 
des sont installees. On admet que la vapeur de chauffage n'est pas 
surchauffee et que le sous-refroidissement du condensat peut etre 
neglige. Le pincement de I'echangeur de chauffage du bouilleur est 
de A 8 p j n B o = 8 K. La temperature de sortie de la solution pauvre du 
bouilleur est alors dg= 125 °C. Par ailleurs, la vapeur de frigorigene 
sortant du dephlegmateur peut etre consideree comme etant dans 
les conditions de saturation. Si nous supposons, la aussi, que la 
chute de pression entre generateur et condenseur est negligeable, la 
pression qui regne dans le generateur est p E = 11,67 bar et 
e 2 = 30”C. 

Evaluation des imperfections de I'absorbeur et du bouilleur 
(generateur) 

Absorbeur : on admettra un ecart des temperatures d'equilibre a 
I'absorbeur sous la pression p F = 2,36 bar de Ad a = 5 K (figure 2). 
Tout se passe done comme si, a la sortie de I'absorbeur, la tempera- 
ture d'equilibre de la solution riche avec I'ammoniac sous la pres- 
sion p F etait 8 a ~ = 8 a + A8 A = 28 + 5 = 33 °C . 

Bouilleur (generateur) : la aussi, on estimera que I'imperfection 
de desorption de la vapeur hors de la solution chauffee sous la pres- 
sion du condenseur pq= 11,67 bar se traduit par un ecart de tempe- 
rature d'equilibre de A 8 q = 5 K (figure 2). Tout se passe comme si la 
solution sortante, bien qu'a la temperature reelle de f? d =125"C, 
n'avait « desorbe » qu'en partie son frigorigene, comme si, ideale- 
ment, elle n'avait ete porte qu'a ( 8 g- A 8 q) = 125 - 5 = 120 °C. 


2.2 Recherche des grandeurs 
thermodynamiques et trace 
des cycles 

Grandeurs thermodynamiques de I'ammoniac pur du circuit CF 

Elies sont donnees dans le tableau suivant. 
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Point 

Pression 

Temperature 

Enthalpie 

massique 

[1] 

Volume 

massique 

[1] 


(bar) 

CC) 

(kJ/ kg) 

(m 3 /kg) 

1 

p F = 2,36 

-5 

1 691,6 

0,532 4 

2 

PC= H.67 

30 

1 709,3 

0,110 7 

3 

PC= H.67 

26 

544,0 

1,663 ■ 10~ 3 

4 

PC= H.67 

26 

544,0 

1,663 ■ 10~ 3 

5 

p F = 2,36 

- 15 

544,0 

(liq. + vap.) 

6 

p F = 2,36 

-11 

1 677,0 

0,517 7 


Nota : la valeur conventionnelle de reference de I'enthalpie massique est de 200,0 kJ/kg 
pour le liquide saturant a - 50 °C. 

La detente 4 — » 5 est isenthalpique. 


Titre en ammoniac de la solution riche x, 

Rappelons : 

— temperature d'equilibre : 33 °C 

— pression correspondant a I'ammoniac a 30 "C : 2,36 bar 

— titre x r [3] : 0,395 

Titre en ammoniac de la solution pauvre x p 

Rappelons : 

— temperature d'equilibre : 120 °C 

— pression correspondant a I'ammoniac a - 15 °C : 11,67 bar 

— titre x p [3] : 0,240 


Grandeurs thermodynamiques des solutions ammoniacales du 
circuit CS 

Elies sont donnees dans le tableau suivant. 


Point 

Titre massique 
en ammoniac 

- . . Enthalpie 

Temperature , on 

r massique 13) 


(kg de NH 3 /kg 
de solution) 

(°C) (kJ/kg de solution) 

a 

x r = 0,395 

28 -104 

b 

x r - 0,395 

? (pour le moment) ? (pour le moment) 
peu differente de 6 a 

c 

x r = 0,395 

~ 88 (point de bulle 175 

sous pq= 1 1,67 bar) 

d 

x p = 0,240 

125 385 

e 

x p = 0,240 

? (pour le moment) ? (pour le moment) 

f 

x p = 0,240 

? (pour le moment) ? (pour le moment) 
mais hf= h e 


Attention ! La valeur conventionnelle de reference de I'enthal- 
pie massique utilisee dans ce tableau est differente de celle du 
tableau de I'ammoniac pur du circuit CF. E lie est ici 0,0 kJ/kg 
pour le liquide a 0 °C. 


2.3 Calculs 


Debit massique de frigorigene (ammoniac), m F en circulation 

Dans I'evaporateur EV, selon la relation (3) comme w= 0 (echan- 
geur statique), on a : 

0 F = m F (h 6 - h 5 ) = m F (h 6 - h 4 ) (kW) (36) 

d'oii : 

0 F 177 

m F = — = = 0,156 kg/s (37) 

(h 6 -h 4 ) (1 677-544) 

Puissance thermique & CD echangee au condenseur 

Le meme debit traverse, en regime permanent, le condenseur. 
D'apres la relation (6), la puissance thermique echangee dans cet 
organe est alors : 

0 C = m F (h 2 -h 3 ) = 0,156x(1 709,3-544,0) = 181,8 kW (38) 

Debit volumique de frigorigene vapeur (ammoniac), V F v, r 
entrant, en 1, dans i'absorbeur 

V FV1 = m F v FVil = 0,156x0,532 4 = 0,083 1 m 3 /s (39) 

Debit volumique de frigorigene vapeur (ammoniac), V F v ,2 
sortant, en 2, du generateur 

V FV: 2 = m F v FV 2 = 0,156x0,110 7 = 0,017 27 m 3 /s (40) 

Debit volumique de liquide frigorigene (ammoniac), V F l, 4 
entrant, en 4, dans le detendeur 

V FL 4 = m F v FL4 = 0,156x 1,663' 10- 3 = 0,259 4- 10- 3 m 3 /s (41) 
Taux de circulation et debits des solutions 

Taux de circulation de la solution riche : d'apres I'equation (10), on 
a, en prenant pour le frigorigene ammoniac qui circule le titre mas- 
sique x FR de 0,995 kg de NH 3 /kg de solution : 

X FR ~ x p 0,995-0,240 

T r r = = = 4,871 

' x r~ x P 0,395-0,240 
kg/s de solution riche par kg/s d'ammoniac (42) 

Debit massique msR de solution riche : 

msR = T c r m F = 4,871 x 0,156 = 0,760 kg/s de solution riche (43) 

Debit volumique Vsr, a tie solution riche sortant de I'absorbeur : 

d'apres la figure 8, la masse volumique de la solution riche a la sor- 
tie de I'absorbeur (x r = 0,395 et 8 a =28°C) est, a peu pres 
Psr ,a = 0,86 ■ 10 3 kg/m 3 . Le volume massique en a est done 
v SR,a = 1/p SR, a ~ 1.163 ■ 10~ 3 m 3 /kg. Le debit volumique de solution 
riche quittant I'absorbeur est done : 

V S R, a = m SR v SR a = 0,760x 1,163 10- 3 = 0,884-10- 3 m 3 /s (44) 

Taux de circulation de la solution pauvre : d'apres la relation (13), 
on a : 

X F R ~ x r 0,995-0,395 

T o, n = = 

x r - x p 0,395-0,240 

= 3,871 kg/s de solution pauvre par kg/s d'ammoniac (45) 
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Figure 8 - Diagramme masse volumique - titre massique 
en ammoniac (source Kaltemaschinenregeln) 


Debit massique msp de solution pauvre : 

msp = t c p mp = 3,871 x 0,156 = 0,604 kg/s de solution 

pauvre ( 46 ) 

Debit volumique Vsp.d de solution pauvre sortant du bouilleur : 

d'apres la figure 8 , la masse volumique de la solution pauvre a la 
sortie du bouilleur (x p = 0,240 et 0 d =125°C) est, a peu pres, 
0,834 ■ 10 3 kg/m 3 . Le volume massique en d est done 
vsp : d = 1,199' 10~ 3 m 3 /kg. Le debit volumique de solution pauvre 
quittant le bouilleur est done : 

V S p,d = m S pV S p d = 0,604x1,199 10- 3 = 0,724 10- 3 m 3 /s ( 47 ) 
Puissance de la pompe PS de solution riche 

Cette pompe deplace le debit volumique de solution riche, Vsr, a , 
de I'absorbeur, ou regne la pression p F = 2,36 bar = 236 000 Pa, vers 
le generateur oil la pression est p E = 11,67 bar = 1 167 000 Pa. Si son 
fonctionnement etait ideal, la puissance mecanique qu’elle consom- 
merait serait (en negligeant les chutes de pression dues a cette 
circulation) : 

Pp S ; . = V S R,a(Pc- Pf) = 0,884 10~ 3 X (1 167-236) 10 3 = 823 W 

( 48 ) 

Mais nous devons tenir compte de I'imperfection du groupe 
motopompe en admettant, pour celui-ci, un rendement energetique 


raisonnable, par exemple Pmp = 0/7. On a alors une puissance reelle 
de cette pompe : 

Vsr, s (P c -Pf ) 0,884- 10- 3 x(1 167-236)- 10 3 

PS,r 9 MP 0,7 

= 1 176 W (49) 

Determination de la temperature de la solution riche a la sortie de 
la pompe (point b) 

En admettant que les echanges thermiques entre la pompe et son 
environnement soient negligeables, alors la variation d'enthalpie 
est, d'apres la relation (18) : 

P PS,r 1,176 

( h b - h a ) = w = — = = 1,547 kJ/kg de solution riche(50) 

msR 0,760 

variation d'enthalpie massique tres faible. E lie entre dans le 
domaine d'incertitude de lecture des enthalpies sur le diagramme 
h-x. D'apres la figure 9, la capacite thermique massique de la solu- 
tion riche dans ces conditions est c ( = 4,47 kJ/(kg • K) . La variation 
de temperature de la solution riche est done 
(9 b -8 g ) = h b -h a / C( = 1,547/4,47 = 0,35 K . Elle est, autant dire, 
negligeable. Sur le trace du cycle des solutions, les points a et b sont 
pratiquement confondus. 

Dans la relation (50), on a utilise la puissance reelle, ce qui 
sous-entend que les pertes de la motopompe sont integrale- 
ment absorbees par le fluide qui circule. Ce n'est pas exact. En 
effet, les pertes du moteur electriques se retrouvent dans I'air et 
non dans la solution en circulation. La valeur trouvee en (50) est 
done une valeur par exces dont I'importance est d'ailleurs 
minime. 

Puissance thermique echangee dans I'echangeur thermique El, 
intersolutions 

D'apres la relation (20), en admettant que I'echangeur El soit par- 
faitement isole thermiquement de I'exterieur : 

0 EI = m SR (h c -h b ) = 0,760 x [175 - (- 1 02,5)] = 211 kW (51) 

(I'enthalpie du point b est prise, pour la forme egale a celle du point 
a = -104 kJ/kg a laquelle on a ajoute 1,5 kJ/kg du a la pompe). 

Refroidissement de la solution pauvre dans I'echangeur El, 
intersolutions 

Les echanges thermiques exterieurs de I'echangeur etant consi- 
deres comme nuls, on a, en appliquant les relations (21) et (22) : 

<P E i = m S p (h d - h e ) (52) 

d'ou : 

h e = h d~— = 385 - _ 35 7 kJ / kg (53) 

m S p 0,604 

Le diagramme h - Ig p nous indique que la temperature de la solu- 
tion pauvre au point e est de I'ordre de 48 ou 49 "C. L'echangeur El a 
done un ecart de temperature de (125 - 88) = 37 K cote chaud et de 
(48 - 28) = 20 K cote froid. Cet echangeur est done parfaitement rea- 
lisable, dans ces conditions. 

On voit I'interet qu'il y aurait a faire un retour arriere de la solution 
pauvre quittant le bouilleur pour prechauffer les plateaux d'epuise- 
ment de la colonne de rectification (cf. figure 6). L'ecart cote chaud 
pourrait ainsi etre diminue, ce qui reduirait l'ecart de temperature 
cote froid. L'absorbeur pourrait etre refroidi a moindre cout. 

Appliquons la relation (2) au regleur de solution pauvre RS. Pour 
ce composant, w= 0 (e'est un robinet, parfaitement statique) et q= 0 


T> 

0) 

o 
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Titre massique en ammoniac x (kg/kg) 

Figure 9 - Diagramme capacite thermique massique - titre massique 
pour I'ammoniac (source Handbuch der Kaltetechnik) 


namique applique a cet ensemble, pris comme circuit ferme, est 
bien verifie. Pour cela, comparons les puissances regues de I'exte- 
rieur et les puissances cedees a I'exterieur, celles-ci devant etre 
egales : 


Puissance regue de I'exterieur Puissance cedee a I'exterieur 
par le composant par le composant 


(kW) 


(kW) 


Evaporateur 

177,0 

Condenseur 

181,8 

Generateur 
[total BO(+) DE(— )] 

300,1 

Absorbeur 

298,4 

Pompe 

1,2 



Total 

478,3 

Total 

480,2 


Le premier principe est done bien verifie a 1,9/479 = 0,004 pres et 
cela malgre I'utilisation d'un diagramme de petite taille. 

La figure 10 montre le trace du cycle suivi dans le diagramme 
enthalpie concentration h-x. 


(il n'y a pas d'echange thermique, la surface d'echange thermique 
offerte etant negligeable). De ce fait, ( hf - h e ) = 0, la detente de la 
solution est isenthalpique et hf = h e . 

Bilans energetiques de I'absorbeur et du generateur 


Note liminaire : pour ecrire les bilans de I'absorbeur et du 
generateur, il nous faut prendre en compte les ecarts des 
valeurs de reference, pour les enthalpies, entre la table de 
I'ammoniac (preferable pour sa precision) et le diagramme 
h-lg p. Pour I'ammoniac pur, la convention tableau (§ 2.2) : 
200 kJ/kg pour I'ammoniac liquide saturant a - 50 "C, donne, 
pour I'ammoniac liquide saturant a 0 "C la valeur de 423,76, 
disons 423,8 kJ/kg alors que la reference du diagramme h-lg p 
est 0 kJ/kg pour I'ammoniac liquide saturant a 0 °C. Nous 
devons tenir compte de cet ecart dans les calculs du bilan en 
diminuant de cette valeur les enthalpies donnees par le tableau 
de I'ammoniac. Ainsi, avec les memes references d'enthalpie 
que pour les solutions (references du diagramme h-x), on a : 

h 1 = 1 691,6 - 423,8 = 1 267,8 kJ/kg 


h 2 = 1 709,3 - 423,8 = 1 285,5 kJ/kg 


Lf) 

CD 

Bilan energetique de I'absorbeur 


o 

CO 

07 

+ mphj- 0 AB - m SR h g + msph f = 0 

(54) 

To 

0 ab = m F h 1 + m S ph f -m SR h a 

(55) 


<P AB = (0,1 56 X 1 267,8)+ (0,604 x 35,7) - [0,760 x (-104)] = 298,4 kW 


Si Bilan energetique du generateur 

O 

o 

- m F h 2 + O m - <P DE - m S ph d + m SR h c = 0 (56) 

0 ) 

~o 

0 ) ... 

( 0 M ~ °de'> = m Fh 2 + m sph d - m SR h c (57) 

o 

® ( 0 m~ 0 DE ) = (0,156 x 1 285,5) + (0,604 x 385) - (0,760 x 175) 

5 = 300,1 kW 

o 

Q. 

£ L'ensemble des calculs des diverses puissances mises en jeu dans 
= les composants du systeme a absorption considere est maintenant 
J termine. II nous faut verifier que le premier principe de la thermody- 

-aj 
CO 

c 
0) 

E 

3 
O 

o 

T> 


Determination des puissances thermiques <t> M et & DE - Epure 
d 'etude de la colonne de rectification (figure 11) 

Sur le diagramme enthalpie-concentration en ammoniac h-x, on 
selectionne les trois isobares qui vont servir a la construction de 
I'epure et qui correspondent toutes les trois a la pression qui regne 
dans le generateur de vapeur d'ammoniac et dans le condenseur : 
Pq = 11,67 bar. Celles-ci sont, rappelons-le : la courbe de rosee R, la 
courbe auxiliaire A et la courbe de bulle B. On porte egalement sur 
le diagramme les trois isotitriques specifiques : celle de la vapeur 
ammoniacale produite x F = 0,995, celle de la solution liquide riche 
x r = 0,395 et celle de la solution liquide pauvre x p = 0,240. 

L' intersection de I'isotitrique x r et de la courbe de bulle B donne le 
point c. La construction classique utilisant la courbe auxiliaire per- 
met de tracer I'isotherme 9 C qui, prolongee jusqu'aux isotitriques x F 
et Xp permet de connaitre les positions limites des deux poles P/im. 
pour la purification et P Em pour I'epuisement. Au point P;, m corres- 
pond, selon le diagramme, I'extraction de chaleur minimale egale a 
Q DF ,jj m . = 252 kJ/kg de frigorigene. On sait qu'il faut adopter une 
valeur plus elevee pour que la droite des poles ait une pente supe- 
rieure a celle de I'isotherme 9 C . On a pris ici, a titre d'exemple, 
Qde= 1,44 x OpEjim . = 252 x 1,44 = 362,9 kJ/kg de frigorigene. On 
obtient ainsi le point P et on peut tracer la droite des poles PP’ qui 
passe par le point c correspondant a la solution riche entrant dans le 
generateur. La longueur du segment P'd donne la quantite de cha- 
leur O m a fournir au bouilleur par unite de masse de solution pauvre 
quittant ce bouilleur. On trouve ici: O m = 592,6 kJ/kg de solution 
pauvre. On peut maintenant calculer les puissances thermiques & M 
a fournir au bouilleur du generateur et & DE a extraire du 
dephlegmateur : 

<P M = q M m S p = 592,4x0,604 = 357,8 kW (58) 

et : 

&OE = Q DE m F = 362,9x0,156 = 56,6 kW (59) 

L'apport energetique global au generateur de vapeur ammonia- 
cale est done : 

0 M~ 0 DE= 357,8 -56,6 = 301 ,2 kW (60) 

Cette valeur est bien conforme a celle trouvee plus haut 
[relation (57), aux inevitables ecarts pres dus aux lectures sur le dia- 
gramme. 
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Figure 10 - Cycle des solutions dans le diagramme enthalpie- 
concentration du couple ammoniac-eau 



Figure 11 - Application de la methode de Ponchon-Savarit au couple 
ammoniac-eau 


Ces resultats etant obtenus, on effectue I'epure de determination 
de la colonne de rectification en suivant les indications donnees au 
paragraphe 1.2. On obtient ainsi I'epure representee figure 11 . La 
colonne de rectification est des plus simples puisqu'elle comporte 
3 plateaux (theoriques), un dans la partie « epuisement », done au- 
dessous de I'arrivee de solution riche, le plateau ©, et deux dans la 
partie « purification », done au-dessus de cette arrivee, les 
plateaux © et ©. 

Coefficient de performance frigorifique de ce systeme 

Le coefficient de performance est : 


Coefficient de performance frigorifique du systeme tritherme 
ideal - Rendement de cycle 

On a: 0 F =-15"C 7> = 258,15 K 

9 C = 30 °C T c = 303,15 K 

8 m = 125 "C T m = 398, 15 K 

Pour un systeme tritherme ideal, on applique [BE 9 734, 
relation (1 )] : 


COP 


a 

F , (id) 


T f (T m - t c ) 

(T c -t f ) t m 


( 62 ) 


COP 


Effet utile 
Ce que I'on paie 


0 P 


177 

357,8+ 1,2 


= 0,493 ( 61 ) 


COP?,, # 


258,15 (398,15-303,15) 

(303,15-258,15) 398,15 


1,37 


Notons au passage que I'energie mecanique que I'on apporte au 
systeme pour le fonctionnement de la pompe est extremement 
faible, pratiquement negligeable devant I'apport thermique au 
bouilleur. 


Le rendement de cycle rj c de notre systeme : 


9c 


COP 

COP° (id) 


0,493 

= 0,36 

1,37 


( 63 ) 
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Valeur assez mediocre, convenons-en. Une recherche d'optimali- 
sation s'impose. 

Dimensionnement des principales conduites de la machine 

Le tableau ci-apres donne les diametres des diverse conduites 
determines a partir des debits volumiques calcules plus haut et des 
vitesses choisies : 



Vitesse 

Debit 

Diametre 

Conduites 

admise 

volumique 

interieur 


(m/s) 

(m 3 /s) 

(cm) 

de frigorigene ammoniac 


Entree de vapeur 
dans I'absorbeur (point 1) 

10 

0,083 1 

10,3 

Entree de vapeur 

dans le condensateur (point 2) 

12 

0,017 27 

4,3 

Conduite de liquide 
(entre les points 3 et 4) 

0,5 

0,259 4 ■ 10~ 3 

2,6 

de solutions ammoniacales : 


Solution riche sortant 
de I'absorbeur (point a) 

1,0 

0,884 ■ 10~ 3 

3,4 

Solution pauvre sortant 
du bouilleur (point d) 

1,0 

0,724 ■ 10~ 3 

3,0 


Pour I'absorbeur : 

0 AB = Veau,ABp gau C gau (8 sAB -9 eAB ) ( 66 ) 

y = 0 AB = 298 ' 4 

eau,AB PsauC eau (e S:AB -e S:AB ) 1 000x4,18x4 

= 0,017 85 m 3 /s soit 64,25 m 3 /h. ( 67 ) 

Pour le dephlegmateur : admettons qu'a I'instar du condenseur le 
rechauffement de I'eau soit de 5 K. On a alors : 


0 DE ~ ^eau, DE P eau c eau(®s, DE~ ®e,DE) ( 68 ) 


Peau c eau(®s, DE~ ®e, DE> 1 000x4,18x5 

= 0,002 708 m 3 /s soit 9,75 m 3 /h. ( 69 ) 

Debit du flulde de chauffage du bouilleur (vapeur d'eau) 

L’enthalpie de condensation de la vapeur d'eau sous 3 bar 
(= 133 °C) est L ve = 2 163,7 kJ/kg [4]. Pour fournir les 357,8 kW de 
chauffage du bouilleur, il faut consommer un debit de vapeur d'eau 
de : 


0 M 357,8 

mvE = = = 0,165 kg/s de vapeur 

LyE 2 163,7 

soit 595 kg/h de vapeur. 


On choisirait les tubes acier, etires sans soudure, dont les dia- 
metres sont les plus proches des valeurs calculees. 

Debits d'eau de refroidissement du condenseur, de I'absorbeur et 
du dephlegmateur 

La capacite calorifique massique de I'eau etant c eau = 4,18kJ/ 
(kg ■ K) et sa masse volumique p eau = 1 000 kg/m 3 , on a les debits 
volumiques suivants : 

Pour le condenseur : 

0 C = V eau , CD P e au c eau( @s, CD ~ ®a, CD ) (04) 

& r 181,8 

Vaail CD — = 

PsauO eau (e s ,CD-e e ,cD ) 1 000x4,18x5 

= 0,008 70 m 3 /s ou 31,3 m 3 /h. (65) 


Remarque finale 

La necessite d'epurer la vapeur ammoniacale pour lui extraire 
la presque totalite de I'eau qu'elle contient impose une colonne 
de rectification. La production du liquide de reflux de cette 
colonne, par le dephlegmateur ou le condenseur dont on 
detourne une partie du condensat, greve le budget energetique 
de la machine de pres de 60 kW (soit environ 16 % de la consom- 
mation d'energie thermique). C'est evidemment un inconve- 
nient du couple ammoniac-eau. 

Malgre I'epuration, un peu d'eau circule dans le circuit frigori- 
fique (0,005 x 0,156 = 7,8 ■ 10~ 4 kg/s ou 2,81 kg/h. On doit assurer 
le retour de cette eau vers I'absorbeur. 

Enfin, a titre comparatif, precisons qu'une machine a com- 
pression monoetagee, de meme puissance frigorifique, 177 kW, 
avec un rendement de cycle de 0,5 et des conditions identiques 
a I'evaporateur et au condenseur, consommerait une puissance 
mecanique d'environ 62 kW et cederait au milieu exterieur la 
puissance thermique de 177 + 62 = 239 kW alors que notre sys- 
teme a absorption en cede (181,8 + 298,4 + 56,6) = 536,8 kW soit 
un peu plus du double. 
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